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ANALOG COMPUTERS 
Laboratory Experiment No. 5 


Solution of Bessel's Differential Equation 


By programming the Analog Computer to solve Bessel's Differential 
Equation, plot J o®) and J, (x) using the Variplotter. Determine from 
the plot the first four zeros of J) and J, (x) and compare with 
tabulated values of these quantities. Note that there is a relationship 
between Jo{x) and J,(x) which should be taken advantage of in this 
problem. 


Discuss the solvtion and any unusual features of your solution. 


Figure 1: Classic Exercise 


Generating Bessel functions 


In June 2022, fellow analog computer enthusiast Dr. CHRIS GILES sent me a 
classic exercise shown in figure 1. This exercise is the basis of this application note. 

BESSEL functions were first described by DANIEL BERNOULLI and later gener- 
alised by FRIEDRICH BESSEL.? BESSEL functions of the first kind are usually denoted 
by J,,(t) and are solutions of the BESSEL differential equation 


Py +ty+  —n7)y=0. (1) 


Sometimes these are called cylindrical harmonics. The parameter n in the equation 
above defines the order. In the following, n = 0 and n = 1 are assumed. 
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For n = 0 equation (1) can be written as 


: 1, 

y= ~ 4 —Y 
after dividing by t? and solving for i. This can be readily transformed into an analog 
computer program by applying KELVIN’s feedback technique. The only thing to take 
into account is the term ¢ which is not well suited for an analog computer due to the 
pole at t = 0. It is far more easy to directly generate the quotient ut as y — 0 with 
1] 0. 

The resulting program is shown in figure 2. Time t has been substituted by machine 
time 7 which is generated using an integrator. The parameter 7 determines how fast 
T rises and should be set so that 0 < 7 < 1 during one computer run. 

The relationship between Jo(t) and J;(t) mentioned in the original exercise is an 
interesting one and can be found in [BRONSTEIN et al. 1989, p. 442] or any other 
standard textbook. In general 


d —n = —n 
S(t In(t)) = OSs lt 


holds, which implies 

I(t) = —Jo(t) 
for the case n = 0. Accordingly J;(7) is readily available in the program as it is just 
~y. 
Figure 3 shows the overall program setup on THE ANALOG THING: A typical 
result is shown in figure 4. Note that the program has not been properly scaled. 7 
was set according to the operate-time of the machine running in repetitive mode. The 
central scaling factor \ was set to get the desired result. © 


3See http: //the-analog-thing. org. 
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Jo(T) 


Ji(r) 


Figure 2: Analog computer program for generating BESSEL functions Jo(7) and Ji(T) 
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Figure 3: Setup for generating Jo(7) and Ji (rT) 
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Figure 4: Typical output for Jo(7) and J;(7) 
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